Slow light, open cavity formation, and large longitudinal electric field on slab 
waveguide made of indefinite permittivity materials 
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The optical properties of slab waveguides made of indefinite permittivity (s) materials (IEM) are 
considered. In this medium the transverse permittivity is negative while the longitudinal permit- 
tivity is positive. At any given frequency the waveguide supports an infinite number of transverse 
magnetic (TM) eigenmodes. For a slab waveguide with a fixed thickness, at most only one TM 
mode is forward-wave. The rest of them are backward waves which can have very large phase index. 
At a critical thickness, the waveguide supports degenerate forward- and backward-wave modes with 
zero group velocity. Above the critical thickness, the waveguide supports complex-conjugate decay 
modes instead of propagating modes. The presence of loss in IEMs will lift the degeneracy, result- 
ing in modes with finite group velocity. Feasible realization is proposed. The performance of IEM 
waveguide is analyzed and possible applications are discussed which are supported by numerical 
calculations. These slab waveguides can be used to make optical delay lines in optical buffers to 
slow down and trap light, to form open cavities, to generate strong longitudinal electric fields, and 
as phase shifters in optical integrated circuits. 

PACS numbers: 42.79.Gn, 41.20.Jb, 04.20.Jb, 78.67.-n 



I. INTRODUCTION 

There is a strong interest in slow light structures and 
devices [1-14]. This is driven by the need for optical 
buffers and optical memory in optical integrated circuits 
and other applications. In order to have small group ve- 
locity, various methods have been exploited to engineer 
material and structural dispersions, such as electromag- 
netic induced transparency [15], coupled resonant struc- 
tures [5, 16], and optical nonlinearity [3]. 

Tsakmakidis et al [6] proposed a novel scheme to real- 
ize trapped rainbow by using negative-index metamate- 
rials. This slow light waveguide requires its core layer to 
be made of double-negative metamaterial (DNM) whose 
permittivity and permeability are both negative. It is 
quite a challenge to realize DNMs with very low loss 
[17] since even moderate loss will destroy the zero group 
velocity mode [18, 19]. Strictly speaking, when loss is 
present as is inevitable in passive systems, storing light 
indefinitely is impossible. However, gain may be intro- 
duced to compensate loss and makes zero group velocity 
possible [12]. Slow light waveguide with single-negative 
materials has also been proposed [12, 13]. 

In this paper we consider planar waveguide made of the 
so-called indefinite medium [20-23] whose permittivity 
and/or the permeability tensors are indefinite matrices. 
For an indefinite medium, the dispersion is hyperbolic 
for one polarization and elliptical for the other. Nega- 
tive refraction, superlens imaging, and hyperlens focus- 
ing [24, 25] can be realized by using indefinite permittiv- 
ity (e) materials (IEMs). Nanowire waveguide made of 



IEM has been considered by Huang et al [9] . These IEM 
waveguides can support both forward- and backward- 
wave modes. High phase index can be obtained for these 
guided modes. Here we further exploit the unique proper- 
ties of IEM waveguide by considering a more manageable 
planar waveguide geometry. We will reveal that these 
waveguides can also support degenerate modes which can 
be used to slow down and trap light, form open cavi- 
ties and generate strong longitudinal electric field. These 
unique properties can lead to a broad range of applica- 
tions. 

The paper is organized as follows. In Sec. II, we 
present the solutions for the TM modes supported by an 
anisotropic slab waveguide. The realization of IEMs will 
be proposed and light coupling into IEM waveguides will 
be discussed in Sec. III. In Sec. IV, various unique fea- 
tures such as slow light, open cavity formation, and large 
longitudinal electric field will be revealed. We conclude 
in Sec. V. 



II. WAVE PROPAGATION IN AN 
ANISOTROPIC PLANAR WAVEGUIDE 

A. Eigen modes equation 

We consider a slab of planar waveguide made of 
anisotropic metamaterial in air. The wave propagation is 
along the ^-direction with phase e % ^ z ~ ut ^ and the trans- 
verse direction is in the x-direction. We consider the case 
of an IEM with 



e z > 0, e x < 0. 



(1) 
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For the transverse magnetic (TM) modes, the magnetic 
field is in the y-direction. The transverse component k x 



2 



of the wave vector inside the metamaterial is 



k x = y/elJ fJ-yk 2 - P 2 /e a 



(2) 



Here kg is the wave number in the vacuum. Since we do 
not consider magnetic materials, we set fi v = 1. Depend- 
ing on the value of e y , this extremely anisotropic medium 
can be uniaxial or biaxial [26-28]. However, since e y is 
not involved with the TM modes, its value is not of our 
interest in this paper. Due to the negativity of e xi the 
dispersion Eq. (2) is hyperbolic instead of elliptic, which 
is shown in Fig. 1. 
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FIG. 1: (Color online) Hyperbolic dispersion relation for the 
FM waves in an IEM. The dashed line (blue) is for a perfect 
electronic conductor (PEC) grating whose transverse compo- 
nent of the wave vector is independent of the longitudinal 
component j3. 



In order to facilitate the analysis of the eigen modes 
and subsequently the calculation of the group index n g , 
we introduce a parameter £ with 



k x = k ^e z (l~l/e x )/^/l-e, 



(G) 



For real e x and e z , one has < £ < 1 for the guided 
modes. For complex permittivity, £ takes complex values. 
We further introduce two parameters 



-l/e a 



(7) 



Note that one has e x = —1/a and e z — cr/g 2 . Thus g 
and a can also be used to characterize the IEM. 
In term of £, a, g, the eigen equation (5) becomes 

^T z k d = F(f, B ,a) (8) 

with the function F defined as 



F 



The phase index is 



2arctan— 
Q 



(9) 



(10) 



Since the planar waveguide is symmetric, the magnetic 
field is 

H y {x) = e ^(x+d/2)^ x < _ rf / 2i 

= Ae lk * x + Be~ ikxX , -d/2 <x< d/2, 
= e -«o(»-d/2) j x > d / 2 . (3) 

Here d is the slab thickness and = \/ f3 2 — k 2 is the 
decay constant in the transverse direction in air. The 
tangential electric field is 

E z (x) = x <_ d/2: 

= — ^(Ae^-Be-*'), -d/2 < x < d/2, 
= _«2 e -«o(*-«J/2) > x > d /2. (4) 

Kq 

The matching of the tangential electric and magnetic 
fields at boundary x = ±d/2 leads to the following eigen 
equation for the TM m modes 



l / v^v n p _1 \ 

kod = — I nnr + 2 arctan — — I . 

(5) 

Here n p = f3/k is the phase index and m = 0, 1, 2, • • •, is 
the parity index of the guided mode. 



B. Asymptotic solutions 

The phase index n p of the guided modes is given im- 
plicitly in the eigen equation (5) or (8). At a fixed wave- 
length or frequency, e x and e z are constant. For real 
values of e x and e z , one can treat the thickness d as a 
function of n p which is chosen as a free parameter. The 
band structure for TM modes on a slab waveguide with 
e x = —3 and e z = 2 is plotted in Fig. 2. 
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FIG. 2: The phase index n p of the guided TM m modes on a 
free-standing planar waveguide of thickness d with e x — —3 
and e z = 2. 

For complex values of e x and e zi the phase index n p will 
also be complex, which is not so easy to solve. One must 
resort to numerical methods. Nevertheless, asymptotic 
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solutions can be obtained analytically, which will reveal 
some of the unique properties of the guided modes on 
IEM waveguide. 

For the modes near the light line, £ — > 0, one has the 
following expansion 



Thus one has 



1 1 

2 + 3g~ 2 - 



£ - ^g(vk d- run) + ^-g 3 (i]k d- mirf 



2 



1 



m ir g g 6 (rjk d - rrvK)\l2) 



with r/ = VI + = v^V^V^i- The phase index 

for the modes near the light line is obtained as 

n p ~ 1 + 1(1 + a- 1 )^ 2 - JfS-a- 1 !^ 

~ l + i( ? 2 (l + o-- 1 )(r7/jod-TO7r) 2 
8 

1 



32 
1 



(1 + <T 1 )m,TTQ 4: (rikod — m7r) s 



with the coefficient & given by 

& = -l + 3cr- 1 - ( T- 2 + (l + cr- 1 )(^- 



--m 7T £1 ) 



(13) 



(14) 



So the modes near the light line is not much different from 
that on a dielectric waveguide [29] . Thus these modes are 
forward-wave modes. 

Due to the hyperbolic nature of the dispersion Eq. 
(2), IEM waveguide supports modes with k — > oo, thus 
n p — > oo. There is no upper bound for n p , unlike the con- 
ventional dielectric waveguide. For these modes, £ — > 1, 
one has the expansion 



(7 V 



m7r + 2arctan 



and subsequently 
^1 



It is evident that validity of £ 
one obtains 



(l + a)e z (hdf 
2(to7t + 2arctang~ 1 ) 2 

1 requires d - 



(15) 

(16) 
0. Thus 



ak d 



m7r+2arctan- 



1 



g-^il + a^kody 



q (1 + g 2 )(mir + 2arctang v ) 2 

(17) 

From this expression of n p , it is evident that the waveg- 
uide supports infinite number of TM modes at any thick- 
ness since all the modes with different parity index m will 
approach the limit n p — > oo as d — > 0. This leads to the 
divergence of the longitudinal electric field E z (see Eq. 
(4)), similar to that on a metallic nanowirc [30]. Large 
longitudinal electric field will be discussed more in Sec. 
IV. 



C. Energy flow and group velocity 



For waves inside the IEM, one has 



(18) 



If k x > y/e^k , j3 will be real and negative if the imag- 
inary part of the permittivity is ignored. So the waves 
inside the IEM will be left-handed, /3S Z < with the 
z-component of the Poynting vector [29] 



1 



Re(E x H* y ). 



(19) 



However for waves confined in the transverse direction, 
S z is no longer uniform. The energy flow of the guided 
modes is better characterized by a quantity P z which is 
defined as [29] 



1 f°° 

- / E x H*dx. 

1 J-oo 



(20) 



Here P z is actually the power per unit length (unit length 
in the y-direction) . For the TM m modes on the planar 
waveguide, P z consists of two parts, one inside and one 
outside the waveguide, 



Pz = P? + P° 



(21) 



with 



pm _ 



pout 



4 E, 



Kl yd ( F 



1 + (-l) m sinc k x d 



4 f \ yr+^ g 2 + e 

\n p \ 1 + {-l) m cosk x d 
4 k q 



n p \ ad 2g y/l — £ 2 



(22) 



The energy flow in the air P° ut is positive as expected. 
However, due to the negative sign of e x , the energy 
flow inside the waveguide is negative and thus contra- 
directional to that in the air. A guided wave is forward 
(backward) wave only if P z is positive (negative). 

The group velocity for the TM modes on the IEM 
waveguide is given by 



v„ = 



dfcp 

d/3 



c 

n a 



(23) 



with c the speed of light in the vacuum and n g the group 
index defined as 



. , drip 
dk 

From the above expression of n p , one has 
dn r , 



"p_ 

dko 



dn p d£ 



drip da 
da dko 



(24) 



(25) 
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with 



dn p 

~W 

drip 



■(n p + —), 

\ an v J 



e 



2n p ^(i-er (26) 

Taking the derivative on both sides of the eigen equa- 
tion (8) with respect to fco, one has 



Fa fco_d£,\ = d£_ + dp „ .da 
o ' dfc dfco 



fco V 2e z dfc, 
Thus one obtains 



+ d„F— . (27) 
dfc 



dfco 
with 



1 



F( 1 + ^sr)-^-^sr 

fco V 2e z dfc / dfc dfc 



(28) 



2e v/T^C 5 



VT+^ g 2 + e i - £ 



2 ' 



2 gy/T^ 



2(l + a)' 



(29) 



After inserting the expressions (25)- (29) into Eq. (24), 
the final expression for the group index n g can be ob- 
tained. 

The expression for v g is very involved. However, the 
group velocity of the modes near light line or with large 
phase index can be easily evaluated. From the asymp- 
totic solutions we have derived in the previous subsec- 
tion, one can see that the guided modes near the light 
line are fast modes with group velocity close to that in 
the vacuum. For the modes with very large phase index 
rip, the group index is also large, leading to the vanishing 
of the group velocity. Due to the finiteness of the ma- 
terial dispersion, the divergence of n g is related to the 
vanishing of d^F which is determined by the geometry 
of the waveguide. One can also verify that d^F = will 
lead to P z = by using the expressions in Eq. (22). 

Using the above expression of P z , one can determine 
the nature of the modes. For the modes near the light 
line, n p — > 1 and £ — > 0, one has P z n ~ — crrf/4 and 
P° ut -> +oo. Thus P z > and the modes are for- 
ward waves. For the modes with n p — > oo and £ — > 1 
which requires d — > 0, one has P z n ~ — oo and p° nt — > 
g 2 /2ko(g 2 + 1). Thus P z < 0, the modes are backward 
waves. For modes with finite n p and n p > 1, P z is finite. 
From the band structure shown in Fig. 2, one would 
expect that P z will decrease with increasing n p . At cer- 
tain critical point d c or n p . c , P z — 0. Further increasing 
n p , P z will become negative, indicating the mode will be 
backward wave. The critical values of the phase index 
n PiC and thickness d c are determined by a critical value 
£ c which satisfies the condition 



2g(i - a ^ . c c 

—— — — = ttitt + 2arctan — . 

ic(Q 2 +i 2 ) Q 



(30) 



The condition of zero group velocity is only determined 
by the parameter g = l/y/—e x e z . The above condition 
can be rccastcd as 



1 + U^- + — ) f tott + 2arctan— )} \ (31) 
2\£ r . P/V pJ\ 



The dependence of the critical values £ c on g is plotted 
in Fig. 3. Once £ c is determined, the critical values of d c 
and be calculated for a given value of a. 




-1/2 
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FIG. 3: The dependence of the critical value £ c for the zero 
group velocity mode on the IEM parameter q = 1/V— e x £ z 
for the TM modes with different parity m. 



From the plotted example (Fig. 2), a few salient 
features are evident. Unlike a conventional dielec- 
tric slab waveguide which supports only forward-wave 
modes, most modes here are backward- wave modes since 
dn p /dd < 0. Only for modes near the light line where 
n p ~ 1, they are forward- wave modes since dn p /dd > 0. 
As shown in Ref. [9] , one can prove that d^F > will lead 
to drip/dd > 0, subsequently P z > and vice versa. Thus 
one can use the sign of dn p /dd to determine whether a 
mode is a forward- or backward-wave mode. Also un- 
like the DNM waveguide which can support single mode 
[31], our IEM waveguide supports infinite number of TM 
modes. That is, it has a very rich band structure for the 
guided TM modes. This can be an advantage if utilized 
properly. Our dielectric cladded IEM waveguide is also 
different from the perfect electronic conductor (PEC) 
cladded IEM waveguide considered in Ref. [32], where 
the supported modes are all backward waves. The coex- 
istence of forward- and backward-wave modes will lead 
to unique properties which will be revealed in the rest of 
the paper. 

So far we have only considered IEM without loss. Any 
realization of IEMs, which will be discussed in the next 
section, will have loss. The definition of IEM will be with 
respect to the real part of the permittivity [27] . However 
all the derivations and expressions will still be valid in 
the presence of loss in IEMs. 
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III. REALIZATION AND LIGHT COUPLING 
TO INDEFINITE PERMITTIVITY MATERIAL 
WAVEGUIDE 

Nature does allow the existence of single-negative ma- 
terials with homogeneous permittivity or permeability 
being negative, such as noble metals, polaritonic materi- 
als, and ferrites. If the underlying crystal is anisotropic, 
the plasmon frequency can be anisotropic, leading to in- 
definite permittivity. One example is bismuth in the THz 
range [32], though bismuth must be placed in liquid he- 
lium temperatures in order to have very small loss. 

Extremely anisotropic materials can be realized easily 
in a broad range of frequencies by proper design of meta- 
materials. One method is through the engineering of arti- 
ficial electric resonances. Another method is through the 
homogenization [27, 33]. For example for a multilayered 
structure of dielectric e a and metal e m , the effective per- 
mittivities can be obtained by using the effective medium 
theory [22, 33], 

e x = fe m + (1 - f)e a , 



Here / is the filling ratio of the metal. For / > / min = 
£a/{£a — Re e m ), one has Re e x < 0. For noble metals 
such as gold or silver, Re e m < in the near infrared 
and visible frequency range. For polar materials [34], 
such as LiTa03 crystal, between the transverse and lon- 
gitudinal optical-mode frequencies, 26.7 THz and 46.9 
THz, the real part of permittivity is also negative. Thus 
if the periodicity of the multilayered structure is much 
less than the operating wavelength, these negative per- 
mittivity materials can be used to design IEMs in their 
respective frequency ranges. 

We design an IEM at A = 1.55/im. The metamate- 
rial is formed by using alternative layers of silver and 
MgF2. At this wavelength, one has e m — —129 + 3.3z 
[35, 36] and e a = 1.9 [37]. Using Eq. (32) with fill- 
ing ratio / = 3.8%, we have e x = -3.0742 + 0.1254i 
and e z = 1.9762 + 0.00003i. This IEM can be used to 
make waveguides of subwavelength thickness. When loss 
is present, the permittivity components in Eq. (1) are 
replaced by the real parts of the permittivity. All of our 
theoretical formulas in the previous section are still valid 
if complex e x and e z are used, which lead to guided modes 
with complex wave number /3 and phase index n p . So- 
lutions must be sought numerically. If loss is ignored at 
this point, the guided TM modes can be easily calculated 
and are plotted in Fig. 4. 

Due to its subwavelength thickness, coupling light into 
an IEM waveguide is a technical challenge. Currently 
there are 4 ways to couple light to waveguides [38]: butt 
coupling, end-fire coupling, prism coupling, and grating 
coupling. The coupling efficiency depends on the cou- 
pling method and the optical properties of the waveguide. 
The end-fire coupling is a butt coupling with a focal lens. 
Multiple modes will be excited in the IEM waveguide by 
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FIG. 4: (Color online) The phase index n p of the first 4 TM 
modes on a free-standing planar waveguide of thickness d with 
e x = -3.0742 and e z = 1.9762 at A = 1.55/zm. The circle 
marks the location of zero group velocity. Insert is for the 
TMo modes. The dashed line is for d — 235 nm. The critical 
thickness for the TMo modes is d c — 238.5 nm. 

butt coupling. However for a waveguide made of realistic 
metamaterial, loss is unavoidable, thus only one or two 
modes will survive over certain distance and eventually 
only one mode will survive after a certain distance. 

Though an IEM waveguide supports infinite number of 
modes, selective excitation of a single mode is possible. 
Thus one can take the full advantage of the rich band 
structure provided by the IEM waveguide. In order to 
excite a single mode, the prism coupling or grating cou- 
pling should be used. However the phase-match condi- 
tion must be satisfied for maximum energy transfer from 
the light source to the waveguide. 

For a simple illustration, we use the prism coupling to 
excite the TM modes in a slab waveguide made of IEM at 
A = 1.55/im. For simplicity, we first ignore the imaginary 
part of the permittivity. Dissipation will be considered 
later in this section. In the range 1.4 < n p < 2.2 and the 
slab thickness d between 225 nm and 240 nm, only the 
TMo modes will be excited as shown in Fig. 4 (see Fig. 
2 for band structure of similar parameters) . The critical 
thickness, such that the forward TMo and backward TMo 
will be merged into a single mode of zero group velocity, is 
d c = 238.5 nm with n p = 1.734. At the thickness d = 235 
nm, two TMo modes are allowed, with n p = 1.541, 1.981. 
The first one is a forward-wave mode while the second is 
a backward- wave mode. 

We place a silicon prism of refractive index n — 3.518 
next to the IEM waveguide. The air gap between them is 
600 nm. We then shine a Gaussian beam into the silicon 
prism. At an incident angle 25.97° inside the prism, the 
forward-wave mode will be excited while at an incident 
angle 34.27°, the backward- wave mode will be excited, 
which are shown in Fig. 5. 

In the above illustration, loss in IEM is ignored. When 
the loss is present, the eigen equation (8) will only allow 
complex solutions. Every mode acquires an imaginary 
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FIG. 5: (Color online) Gaussian beam excitation through 
prism coupling of the forward-wave (a) and backward-wave 
mode (b) on an IEM waveguide at incident angle 25.97° 
and 34.27°, respectively. The IEM has lossless permittivity 
e x = -3.0742 and e z = 1.9762 at A = 1.55/im. The air gap 
between the prism and the waveguide (d — 235 nm) is 600 
nm. Plotted is the absolute value of the magnetic field H y . 

part. There is no clear cut distinction between propagat- 
ing modes and decay modes. Duo to the presence of loss 
in the IEM, the forward-wave TMq modes will not be 
degenerate with TM^ at any thickness. This is shown in 
Fig. 6 in the complex-^ plane. 




0.2 0.4 0.6 0.8 1 

Re^ 

FIG. 6: (Color online) The splitting of the forward- wave (red) 
and backward-wave (blue) TMo modes in the complex-^ plane 
in the presence of loss in the IEM with e x = -3.0742 + 0. 1254i 
and e z = 1.9762 + 0.00003i at A = 1.55/im. The dot-dashed 
line (green) is the solution for TMo modes on the IEM waveg- 
uide without loss. 

The splitting of solutions in the complex-^ plane will 
lead to the split of the phase index n p which is plotted 
in Fig. 7. Consequently, there is no well-defined critical 
thickness. 

At a thickness d — 235 nm, the two TM modes 
have n p = 1.5190 + 0.0672i, 2.0026 - 0.1437i for the 
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FIG. 7: (Color online) The complex phase index n p for the 
forward-wave (top) and backward-wave (bottom) TMo modes 
on an IEM waveguide with e x = -3.0742 + 0.1254i and e z = 
1.9762 + 0.00003? at A = 1.55pm. The dot-dashed line (green) 
is the phase index of the propagation TMo modes on the IEM 
waveguide without loss. 



forward- and backward-wave modes, respectively. Thus 
the forward- and backward-wave modes will propagate a 
finite distance before lose most of their power. Since the 
forward- wave mode is closer to the light line, thus the loss 
will have a smaller effect while the backward- wave modes 
corresponding deeper penetration into the ohmic metal, 
thus will have shorter decay length. The excitation of 
these two modes is illustrated in Fig. 8. 



IV. APPLICATIONS OF WAVEGUIDE MADE 
OF INDEFINITE PERMITTIVITY MATERIALS 

A. Phase shifter and longitudinal electric field 
enhancement 

One salient feature of the modes on the IEM waveg- 
uide is the large phase index. The high phase index is 
due to the hyperbolic dispersion in the metamaterial. 
Nanowires waveguide based on these slab waveguides can 
be used for phase shifter with small footprint in optical 
integrated circuits. The presence of loss in the metama- 
terials will restrict the use of the waveguide. However for 
many applications other than the long-haul transporta- 
tion, short waveguides have the advantage of small sizes. 

As examples, we consider a slab waveguide made of an 
IEM with e x = -3.0742 + 0.1254i and e z = 1.9762 + 
0.00003i at A = 1.55/im. The realization of this IEM was 
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FIG. 8: (Color online) Same as Fig. 5 except that the IEM has 
loss with e x = -3.0742 + 0.1254i and e z = 1.9762 + 0.00003i. 
The air gap is now 400 nm. 



discussed in the previous section. If no loss is present, the 
waveguide made of this IEM will support forward-wave 
and backward-wave modes as shown in Fig. 4. 

At a thickness d = 160 nm, the waveguide supports 
two TMo modes, one is a forward-wave mode with n p = 
1.1211 + 0.0035i and the other is a backward-wave mode 
with n p = 4.0298 - 0.1200z. The TMf with phase index 
Re rip ~ 4 can be used to create large phase shift over 
a short distance. In order to have a phase shift of it by 
this slab waveguide, the propagating length can be about 
192.3 nm. Over this distance, the damping is about 1.26 
dB. For excitation of such high phase index modes, grat- 
ing coupling should be used. 

Besides the large phase index, modes with large longi- 
tudinal electric field can be excited on the IEM waveg- 
uide. Recently there is a strong interest in large longitu- 
dinal electric fields [39]. Large longitudinal electric field 
can have a lot of applications, such as superfocusing to 
beat the diffraction limit [40, 41] and trapping metallic 
nanoparticles in optical tweezer [42]. 

Due to the confinement in the transverse direction, 
there is a 7r/2 phase difference between E z and E x in 
air. The transverse electric field is 

EJx) = A e «o(x+d/2) j x <_ d / 2; 

fc 



p 



£ 



(Ae ik °> x + Be- ik * x ), -d/2 <x< d/2, 
e -«o(x-d/2) j x > d /2. (33) 



For the waveguide thickness d = 160 nm, the electric field 
of the two TM modes are plotted in Fig. 9. 
Here we consider the ratio 

s = \E z \ max /\E x \ max . (34) 

One can see that s ~ 50% for the TM^" mode and s ~ 
100% for the TM^ mode. The ratio is comparable and 



even stronger than that on a silicon nanowire waveguide 
[39, 43]. If one shrinks the thickness, even larger ratio 
would be expected. The fact that the IEM nanowire 
waveguide will have large longitudinal electric field is due 
to the hyperbolic dispersion Eq. (2), which allows much 
stronger confinement of light in the transverse direction, 
thus stronger longitudinal electric field. 
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FIG. 9: (Color online) The longitudinal and transverse elec- 
tric fields of the two TMo modes on a slab waveguide of 
thickness d = 160 nm with e x = -3.0742 + 0.4896i and 
e z = 1.9762 + 0.00003i at A = 1.55pm. One has n p = 
1.1211 + 0.0035i, 4.0298 - 0.1200i, respectively. 

We further derive analytical expression for the ratio s. 
One has 



\E X 



z p |max(l>v/l + £7e 2 |)- (35) 



and 



\E, 



m = 0, 



m > 0. (36) 



By using the above expressions and also the eigen equa- 
tion (8), the ratio s can be calculated for modes with 
different waveguide thickness d as shown in Fig. 10. De- 
pending on the value of e x and e z , the ratio s can be over 
100%. 



B. Light wheel and open cavity formation 

Resonances are ubiquitous. They have many applica- 
tions such as storing and confining energy, enhance the 
field concentration, and improve the detection accuracy. 
To have a resonance, a compact space is required, such as 
cavities. Examples are the quantum dots, microwave cav- 
ities, photonic crystal microcavities. Recently, negative- 
index metamaterials are also used to form open cavities 
[44-47] , such as the checkerboard open resonators [48] . 
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FIG. 10: (Color online) (a) The real part of the complex 
phase index n v and (b) the ratio s of the longitudinal and 
transverse electric fields of the TM modes on a free-standing 
slab waveguide of thickness d with e x = —3.0742 + 0.4896i 
and s z = 1.9762 + 0.00003i at A = 1.55/um. 



For optical integrated circuits, one of the most diffi- 
cult tasks of nanofabrication is the alignment of different 
parts and devices. It would be more desirable to have a 
open cavity at any location. 

Recently, a new concept, a light wheel, has been de- 
veloped [49]. This is formed in a composite waveguide, 
which is made of an ordinary slab waveguide coupled 
with a properly designed slab waveguide made of DNM. 
If these two waveguides are separated infinitely away, the 
ordinary waveguide support a single forward-wave mode. 
The DNM waveguide supports a single backward-wave 
mode with the same phase index. Once these two waveg- 
uides are placed in the vicinity of each other, the com- 
posite waveguide no longer supports propagating modes 
at the same wavelength. Instead it will support complex- 
conjugate decay modes. These two decay modes will form 
the so-called light wheel [49] . 

In order to have complex-conjugate decay modes, the 
waveguide should first be able to support degenerate 
propagating modes. In the example we considered in Sec. 
Ill, the slab waveguide is made of a metamaterial with 
e x = -3.0742 and e z = 1.9762 at A = 1.55/xm. Below the 
critical thickness d c = 238.5 nm, the waveguide supports 
two modes of different phase indices (see Fig. 4). One 
mode is a forward-wave TM mode and the other is a 
backward-wave TMo mode. At the critical thickness, the 
effective thickness of the waveguide is zero due to the neg- 



ative Goos-Hanchen lateral shift and a double light cone 
will be formed [6]. However above the critical thickness, 
the waveguide supports no propagating modes. Instead, 
it supports complex-conjugate decay modes. For exam- 
ple at d = 239 nm, one has n p = 1.7286 ± 0.0373?. When 
an incident beam with /3 — 1.7286fco hits the waveguide, 
the two decay modes will be excited, one decays along 
/3 and the other in the opposite direction of /?, thus an 
open cavity will be formed. This cavity can be formed at 
any location along the IEM waveguide, which is shown 
in Fig. 11(a). 
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FIG. 11: (Color online) Light wheel formation through a 
prism coupling on an IEM slab waveguide of thickness d — 239 
nm at A = 1.55/i. (a) Lossless IEM with e x ~ —3.0742 and 
e z = 1.9762 and (b) IEM with e x = -3.0742 + 0.4896i and 
£ 2 = 1.9762 + 0.00003i. The incident angle of the Gaussian 
beam in the silicon prism is 29.4°. The air gap between the 
prism and the waveguide is 400 nm. Plotted is the magnetic 
field H v . 

When dissipation is present, the decay modes no longer 
form a complex-conjugate pair. However the imaginary 
part of their complex phase indices will still have oppo- 
site signs, thus light wheel formation is still allowed. In 
our example as shown in Fig. 7, for d = 239 nm, the 
two modes has the complex phase index n p — 1.5877 + 
0.1282z, 1.8730 - 0.2042. The formation of light wheel is 
shown in Fig. 11(b). 



C. Slow light waveguide 

For the anisotropic waveguide we have considered, 
e x < 0, one has P] n < and P° ut > if one sets > or 
n p > 0. If P z = Pf + P° ut < 0, the mode is a backward- 
wave mode since the total energy flow is opposite to the 
phase velocity. Otherwise, the mode is a forward-wave 
mode. Unlike the dielectric waveguide or the PEC grat- 
ing, the IEM waveguides support both forward and back- 
ward waves. At the critical thickness d Cl the backward 
and forward modes become degenerate, the energy flow 
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inside the waveguide cancels out that in the air. One 
can prove that at d c where P z = 0, the group velocity 
is indeed zero. One does not need to know the mate- 
rial dispersion to locate the zero group velocity point. 
This is due to the fact that for these waveguides, the 
dispersion due to geometric confinement dominates the 
material dispersion at and around the critical thickness. 

This unique property of the modes on an IEM waveg- 
uide can be used to slow down and even trap light. Even 
though the waveguide supports infinite number of TM 
modes at any fixed thickness and frequency, with appro- 
priate laser coupling, the excitation of the TM m modes 
with m > 1 in the waveguide can be suppressed or even 
eliminated. Due to the material dissipation, the first TM 
mode will propagate the longest distance. The rest of 
the TM modes will all decay out at about half the de- 
cay length of the first TM mode. It is the TM band 
which can be used for slow light application, though other 
TM modes can also be used. Unlike the double negative 
waveguide [6], the IEM waveguide will slow down and 
trap light if one increases the thickness to the critical 
thickness d c . These waveguides can thus be used as de- 
lay line in optical buffers [5]. 

For the butt-coupling, multiple modes are excited. To 
excite a single mode, prism coupling or grating coupling 
can be used. By tapering the waveguide thickness, any 
branches of the guided modes can be accessed in princi- 
ple. 

For an asymmetric waveguide where an IEM waveg- 
uide is placed on a dielectric substrate, the excitation of 
some unwanted TM modes can be avoided. A sketch of a 
tapered IEM waveguide on a glass substrate is illustrated 
in Fig. 12. Let the IEM have e x = -3.0742, e z = 1.9762 
and the glass substrate have refractive index n — 1.5 at 
A = 1.55 fj,m. The phase index n p as the function of 
thickness d is plotted in Fig. 13. Due to the presence 
of substrate, the phase index of the guided TM modes 
starts at n p — 1.5 with do — 116.5 ran. The critical 
thickness of the TMo modes is d c — 174.4 nm. Thus if 
the initial thickness of a tapered IEM waveguide is less 
than 115.8 nm, the excitation of the forward-wave TMo 
mode can be avoided. If one further ignores the excita- 
tion of TM m with m > 1 due to the loss or the coupling 
strength, the tapered waveguide will allow only a single 
mode in practice. If the thickness is gradually increased 
to d c = 175.6 nm, the wave will be stopped there. 

However, loss is ubiquitous. The presence of loss in the 
IEM will prevent forward-wave modes and the backward- 
wave modes from joining together as shown in Fig. 7, 
leading to finite group velocity. For the IEM we have 
designed in the previous section, one has the material 
dispersion for the IEM 
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FIG. 12: (Color online) A sketch of tapered slow light waveg- 
uide made of IEM on a glass substrate. 




0.5 1 
waveguide thickness d (nm) 

FIG. 13: (Color online) The phase index Tip SbS Sb function of 
the waveguide thickness d of the first 4 guided TM modes 
on an IEM slab waveguide placed on a glass substrate with 
refractive index n = 1.5. The IEM has e x = —3.0742 and 
e z — 1.9762 at A = 1.55 /im. The circle marks the location of 
zero group velocity. Insert is for the TMo modes which starts 
at n„ — 1.5 with d = 115.8 nm. 



dcr 

dfc 



/ de m 
dk 



(37) 



Here we have ignored the dispersion of e a in Eq. (32). 
We further assume a Drude model dispersion for silver. 
Thus one has 



de r . 



dk 



1 

feo 



2(1 - e m ) - ilme, 



1 



1 



(38) 



For the exemplary IEM waveguide we have designed 
in the previous section, the group index of the modes 
is calculated and plotted in Fig. 14 by using the above 
dispersion of the IEM. One can clearly see that due to 
the presence of loss in the permittivity, especially in e x , 
the group index is greatly reduced. In the example we 
considered, we get n g ~ 10. The fact that the group 
index n g for the backward-wave TMq 3 modes is negative 
is due to our convention of adopting positive phase index 
n p throughout the paper. 

However gain can be introduced so that at a critical 
gain, the two branches of modes will recombine at a crit- 
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FIG. 14: (Color online) Group index n g for the forward- 
wave (top) and backward-wave (bottom) TMo modes on 
a free-standing IEM with g = 0.4055 + 0.0083i and a = 
0.3247 + 0.0132i at A = 1.55 fim. The dispersion of the IEM 
is discussed in the text. 

ical thickness with the same real phase index. One can 
coat gain medium such as semiconductor optical ampli- 
fiers (SOAs) [50] by IEM. The gain in the SOA can com- 
pensate the loss in IEM and zero group velocity mode 
can be recovered. Recent progress has also been made 
to use gain to boost the propagation of surface plasmon 
polaritons [51]. 

We point out that the effective medium theory for the 
metal and dielectric multilayered structure is still valid 
when one goes down in frequency from the visible, near 
infrared to THz waves and microwaves. The multilayered 
structure will still be able to support zero group velocity 
modes when one decreases the frequency. But the phase 
index when the group velocity is zero will increase. In 
the limit of perfect electronic conductor, e m — » oo, one 
has 

e x =oo, e z =e a /(l-f). (39) 

In order to give correct dispersion for waves in the meta- 
material, k x = ^fe^ko for the TM modes (see Fig. 1), 
one should have 

to = 1 - /. (40) 

A free standing slab with the above permittivity and per- 
mittivity will support guided TM modes. Waveguide can 
also be formed by sandwiching this metamaterial with di- 
electric on both sides, or dielectric on one side and PEC 
on the other. In the case of a PEC grating with finite 
depth, the above effective indices will give rise to the 



so-call spoof surface plasmons [52] . Only when the phase 
index goes to infinity, zero group velocity can be reached. 
This leads to the fact that the supported modes are all 
forward- wave modes since the energy flow inside the PEC 
grating is zero and that in the air are along the phase ve- 
locity. Though this type of waveguides can support TM 
modes with very small group velocity [11, 53], they can 
not be used to stop light. Even in the THz range, the 
phase index when the group velocity is zero will be very 
high that in practice, a metallic grating structure will not 
be able to stop light. Due to the hight phase index of the 
guided TM modes, strong reflection is expected along a 
tapered metallic grating. Only when the frequency is 
above the collision frequency, which is about 10 THz for 
silver [37] , the dispersion of a metallic grating will be hy- 
perbolic and the grating structure may be used to stop 
light. 



V. CONCLUSIONS 

In this paper, we consider wave propagation in a slab 
waveguide with anisotropic optical constant. For ex- 
tremely anisotropic cylinder where the transverse com- 
ponent of the permittivity is negative and the longitudi- 
nal is positive, the waveguide supports infinite number of 
TM modes. Among the supported TM modes, at most 
only one mode can be forward wave. The rest of them 
are backward waves. 

Possible realization of these IEM waveguides are pro- 
posed by utilizing alternative layers of metal and dielec- 
tric. Light couplings to IEM waveguides have also been 
discussed. To take full advantage of the rich band struc- 
ture provided by the IEM waveguide, prism coupling and 
grating coupling can be used to selectively excite the 
guided modes. 

Four unique properties have been revealed for the 
modes on slab waveguides made of IEMs. The first is 
that the backward- wave modes can have very large phase 
index. These waveguide can be used as phase shifters in 
optics and telecommunications. The second is the large 
longitudinal electric field of the modes due to the hy- 
perbolic dispersion of the metamaterial. The third is 
the formation of open cavity along in the waveguide due 
to its support of complex-conjugate decay modes above 
the critical thickness. The forth is that the waveguide 
supports modes of zero group velocity. This is due to 
the fact that the waveguide can support both forward- 
and backward-wave modes at a fixed thickness. If the 
waveguide is tapered, at certain critical thickness, the two 
modes will be degenerate and carry zero net energy flow. 
At other thickness, these waveguides support modes with 
small group velocity. The presence of loss in IEMs will 
destroy the zero group velocity modes but gain can be 
introduced to compensate loss. These waveguides can 
thus be used as ultra-compact delay line in optical buffers 
[5]. Many more interesting features many further emerge. 
The above features can lead to potential applications of 
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these waveguides in optical integrated circuits. 

In this paper, only indefinite e mctamaterials are con- 
sidered for waveguide applications. Similar results can 
also be obtained for waveguide made of indefinite perme- 
ability (/i) metamaterials (IMMs). As the case for any 
application of metamaterials with negative parameters, 
loss will severely deteriorate the performance. The pro- 
posed design of IEM is by no means optimal. Similar 
to the IEM waveguide, IMM waveguide will have lim- 
ited application due to large loss. A recent design of the 
split-ring resonators used in Ref. [54] for cloaking has 
very small loss of the permeability and can be used to 



realize the applications we have revealed in this work. 
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